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Sjoqvist, Pati, Ekert, Anandan, Ericsson, Oi and Vedral (Phys. Rev. Lett. 85, 2845 [2000]) have 
recently "provided a physical prescription based on interferometry for introducing the total phase of 
a mixed state undergoing unitary evolution, which has been an elusive concept in the past" . They 
note that "Uhlmann was probably the first to address the issue of mixed state holonomy, but as a 
purely mathematical problem" . We investigate possible relationships between these "experimental" 
and "mathematical" approaches, by examining various quantum-theoretic scenarios. We find that 
" the two methodologies, in general, yield inequivalent outcomes. 
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^ ! I. INTRODUCTION 

(N ; 

In the introductory paragraph of their recent interesting letter, Sjoqvist et al state that "Uhlmann was probably 
\ the first to address the issue of mixed state holonomy, but as a purely mathematical (emphasis added) problem. In 
■ contrast, here we provide a new formalism of the geometric phase for mixed states in the experimental (emphasis 
added) context of quantum interferometry" jij. Left unaddressed there, however, was the nature of any possible 
relationship between the two approaches. This question merits particular attention since Dittmann has "shown 
that the connection form (gauge field) related to the generalization of the Berry phase to mixed states proposed by 
Uhlmann satisfies the source-free Yang-Mills equation *D * Duj = 0, where the Hodge operator is taken with respect 
to the Bures metric on the space of finite-dimensional nondegenerate density matrices" 0. This fact, Dittmann and 
Uhlmann noted, "may be seen as an extension to mixed states of numerous examples relating the original Berry phase 
to Dirac monopoles, and the Wilczek and Zee phase, to instantons" ||, sec. 1]. In their work, Sjoqvist et al do express 
O |. the geometric phase in terms of an average connection form |l], eq. (16)] (for which no Yang-Mills characterization has 
been given). Zee has reported that the non-Abelian gauge structure associated with nuclear quadrupole resonance 
does not satisfy the source- free Yang-Mills equation B (cf. ||). In in line with certain work of Figueroa-O'Farrill 
concerning calibrated geometries [Q, sec. 6], we have begun a study of the "Dittmann-Bures/ Yang-Mills" field over 
the eig/it-dimensional convex set of i/iree-level density matrices (making use of the formulas for the Bures metric 
developed in ||, based on the parameterization in ||) by finding certain self-dual four- forms (ip) with respect to the 
■ ■ »h , Bures metric. The eigenspace decomposition of the associated endomorphism ip was found to yield four octets and 
rS \ three singlets as factors 

The research presented below represents our effort to gain a formal understanding of the relations between the 
"mathematical" methodology of Uhlmann and the "experimental" approach of Sjoqvist et al. We hope that this will 
prove useful, among other things, in the devising of a physical scheme for implementing the concepts of Uhlmann (cf. 
@ sec. 4]). 

In the interferometric/experimental case Q, each pure state that diagonalizes the initial density matrix (p(0)) is 
parallel transported separately from other distinct diagonalizing pure states, so one simply weights or (incoherently) 
averages the individual pure state (single-state) interference profiles by the corresponding eigenvalues of p(0) to obtain 
the mixed state results . In the notation of Q, the total geometric phase is the argument of '^2 k Wki / ke l ^ k , and the 
visibility is its absolute value, where Wk, Vk, 4>k are, respectively, the eigenvalue, visibility and geometric phase obtained 
for the fe-th constituent pure state of p(0) fiL eqs. (10), (15)]. The phase introduced by Sjoqvist et al "fulfills two 
central properties that make it a natural generalization of the pure state case: (i) it gives rise to a linear shift of the 
interference oscillations produced by a variable U(l) phase and (ii) it reduces to the Pancharatnam connection for 
pure states" [Q . The mixed state generalization of Pancharatnam's connection — which asserts that p(t) and p(t + dt) 
are in phase if Tr[p(t)U (t + dt) f U(t)} is real and positive — can be met only when |l], eq. (12)] 

Tr[ P (t)U(t)U\t)}=0. (1) 

This is the parallel transport condition of Sjoqvist et al for mixed states undergoing unitary evolution. 
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In the approach of Uhlmann ][Ll) the parallel transport takes place in a larger Hilbert space, defined by purification. 
The geometric phase of Sjoqvist et al can also be understood using purification pj], though the unitary operator that 
acts on the ancilla or auxiliary state there is unconstrained and can in fact be anything, such as the identity operator. 
A "parallel transport of a density operator amounts to a parallel transport of any (emphasis added) of it purifications 
@, P- 2848]. 

In infinitesimal form, the parallelity condition of Uhlmann for a density matrix p(t) reads eq. (19)], 

W(tyW(t) = hermitian, (2) 

where W(t) = p(t)V{t). The unitary operator V(t) can be interpreted as acting on the ancilla in the purification 
of p(t) and is determined by this parallelity condition (^]). For example, if the system, initially in the state p(0), 
evolves unitarily under the evolution operator U(t) one may write W(t) = U(t)y/p(0)V(t). The Uhlmann parallelity 
condition, then, reads 

jmu\t)U{t)JW) = V(t)V(tfp(0) + p(0)V(t)V(t?. (3) 

With given input p(0) and U(t) this equation needs to be solved for V(t). The geometric phase 7 S for any specific T 
can then be computed using the mixed state Pancharatnam connection |jj , 

lg (T) = arg Tr(Ut {T ) y /rffi u i (T) ^q)) . (4) 

In the three sets of analyses reported below (sec. II), we will implement certain formulas of Uhlmann |Tl],[l2|]. We 
compare the two methodologies under consideration in the context of: (a) unitary evolution over geodesic spherical 
triangles, the vertices of which are associated with certain specific mixed two-level systems; (b) unitary evolution of 
such two-level systems over certain circular paths (0(3)-orbits); and (c) unitary evolution over such circular paths of 
fc-level Gibbsian density matrices |T^,[l3| . In case (a) we find an interesting correspondence between the two approaches 

for mixed states having Bloch vectors of length ^/| and in (c) that the results of the Uhlmann methodology converge 

to those of Sjoqvist et al if a certain variable (a), a function (sech^-) of the inverse temperature parameter (a), is driven 
to zero. If instead of strictly mixed states, we consider the unitary evolution of pure states, then both methodologies 
simply reproduce the corresponding Berry phase for each of the three scenarios. However, for general mixed states, 
the two methodologies typically yield inequivalent results. 



II. ANALYSES 



A. Unitary evolution over geodesic spherical triangles 

For a qubit (a spin-i particle), whose density matrix can be written as, 

p(0) = ~(l + rr-<7), (5) 

where f is a unit vector, r is the constant (length of the Bloch vector) for unitary evolution and a is a vector of the 
three (non-identity) Pauli matrices, Sjoqvist et al obtained a formula for the geometric phase jj], eq. (25)], 

7 9 [r] = — arctan^rtan^ . (6) 

Here T is the unitary curve in parameter space traversed by the system and the solid angle subtended by the natural 
extension of T to the unit (Bloch) sphere. This formula reduces to — ^ for pure states (r = 1), as it, of course, should. 



The visibility was expressed as jj], eq. (26)] 



v = rj\J cos 2 — + r' 2 sin 2 — . (7) 

For the case of cyclic evolution, which will be our only area of concern here, rj = 1. 

We focus on the case of cyclic evolution for geodesic triangles, in particular, since Uhlmann has an enabling formula 
"relevant for the intensities", |11, eq. (24)] 
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(8) 



, 2/-, x2°12 + a 23 + 0,31 + 6 , ,3 . x4 

+r (1 — rj 1- 4r(l — r) + (1 — rj . 

Here 

ou := (1|2)(2|1), a 32 i = (3|2)(2|1)(1|3), . . . , (9) 

where gives the projection operator Pj associated with the j-th density matrix (pj) in the geodesic triangle, 

that is, 

p,-®^P&. do) 

Uhlmann also apparently has an explicit formula for geodesic quadrangles, but only the leading terms of it have been 
published [[LlJ eq. (25)], that is, 

V a i2a23a34a 4 iTrp(0)?7 = r 5 a i32 i+ (11) 



+r 4 (l - r) 



6&4321 - 0432 - a 42 i + 3a 43 i + 3a 32 i + 043 + 0-21 + 041 + a 2 3 + a 42 - 3ai 3 + 2 



Without loss of generality, let us set one of the three vertices of the geodesic spherical triangle swept out by the 
unit vector f (||) in the unitary evolution to be (0,0,1). We use the standard spherical coordinate form for a point on 
the sphere, that is (sin$cos0,sin#sin</>, cos#), and take the other two vertices to be parameterized by the angular 
pairs (9i,(f>i) and (^2; ^2): so the preassigned vertex (0,0,1) corresponds simply to 81 = 0. Then, using the notation 
4> = 4>i — 4>2 , we have for the solid angle uM , 



£1 = — 2arccot ( cot <f) + (cot 61 + esc 61) (cot 62 + esc 62 ) esc <p \ . 



(12) 



(This formula, as well as all the succeeding ones here, depend implicitly or explicitly upon <f> and not on <f>i and 4>2 
individually.) Let us also denote 



fi = (1 + cos6*i)(l + cos# 2 ) + cos 4> sin B\ sin 6 2 



(13) 



The parameter /i is confined to the range 4], for < 9%, 62 < 7r, < (f>i, <j>2 < 27r. It is equal to its minimum -~ 
at 0i = #2 = and ([> — tt; equal to its maximum 4 for #1 = #2 = and equal to zero if either Oi — ir or 9 2 = n. 

Making use of (|8|) and (^|), we have for the "Uhlmann geometric phase" for the general geodesic triangle specified 
that 



7 5 = arg(TrpiC/) = arctan^J , 



(14) 



where 



a = r sin <p sin t>i sm O2 , 



/? = 4+( M -10)r 2 +6r 4 . 



(15) 



We found that 



so one can write @ 



tan ■ 



n 



= — arctanfrtan — ^ = arctanf — ^) 
\ 2 J \fj,r 2 J 



7s 



(16) 



(17) 



as well as 
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7 g = arctan^— J — arctan^ — —J. (18) 

Let us note the exact relations between the two forms ("mathematical" and "experimental") of geometric phase, 

tan( 7 J 4-10r 2 + 6r 4 8 - ur 2 (3 

f-^- = l + s =1 + - ^ = — 

tan(7 g 



which are unity for both r = 1 and r = J | . ( Anandan et al Jl5| assert that "a fully gauge- invariant description of 
the geometric phase requires that it be defined modulo 2ir v . Of course, tamp = tan (ip + it).) For r = 1 — Ar, 



tan7 9 / 4Ar 
tan7 ff 



(20) 



For r = w | + Ar, 



tan7 g 

For r — Ar (that is, a neighborhood of the fully mixed [classical] state, r = 0) 



tan 7g _ ^ 2V6Ar \ 



2^^^^1 + i^!:. (22) 
7 5 tan7 9 /iA^r 



The "Uhlmann visibility" is 



/ a 2 + B 2 

v = \Ty Pi U\ = J / 2 , (23) 



which is precisely unity for r = 1, as is also the case for v, given by (^). Now, 



v 1 I a 2 + /i 2 r 4 



v r 2 y a 2 + B 2 



(24) 



This ratio is always unity for r = 1 and always | for r = i / 1. (Recall, on the other hand, that the ratio (19) of the 



tangents of the two geometric phases is unity for r — y f ■) For r = 1 — Ar 

v ' V ^ (oyV 3 ) 2 + ji 2 



- 1 • >Ci • !'' .., )• i- >r ') 



For r = + Ar, 



For r = Ar, 



In the limit r — > 0, 



^ 6 m } y (26) 

zv 2\ v 2(a/r 3 ) 2 + 3^ 2 / 
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B. Unitary evolution over certain circular paths 



Now let us similarly study the unitary evolution of spin-i systems over circular (0(3)-orbits) — as opposed to 
geodesic triangular — paths. We utilize the analyses of Uhlmann in |l2). Again, we start with an initial density 
matrix p(0), the Bloch vector of which is proportional to (0,0,1), while the axis of rotation about which the state 
evolves is it = (0, sin£, cos£). The resulting curve of density matrices is given by 

4> -> U(4>)p(0)U(-<t>), U{<j>) = e-^m^+cose/^ (29) 

1/2 

and the associated parallel lift of this curve with initial value p(0) by 

U{<j>)p{0) 1,2 V{<j>), V{<j>) = e #(^s4+™ S «^) ) (30) 

where 

2 P (0) 1 / 2 J yP (0) 1/2 = P (0)Jy + J yP (0), (31) 

and the J's are, of course, generators of an irreducible representation of SU(2). We have the relation 

e = cos- 1 f-^y (32) 



2tt, 

The associated holonomy invariant for a complete rotation about ~ft is 

p(0) 1/2 C/(27r)^(2^)p(0) 1/2 . (33) 

(In jl3], sec. V], we have studied certain "higher-order" holonomy invariants.) We then have the result (using the 
" " " notation now, rather than " ~ " as before, to refer to results obtained based on the unitary circular evolution 
analysis of Uhlmann Jl^] ) 

tan 7s TrytanKCothTry 

^- = X sec £ coth 7tx tan ( 7r cos £) = , (34) 

tan 7 S k 

where 



r 1 — 2 — r 2 cos 2£ 

o » ( 35 ) 



and 



7rV-l + r 2 -y 2 Q 

— =ttcos4 = - — . (36) 



We plot this ratio (||) in Fig. [| 




FIG. 1. Ratio (Q) of tangents of geometric phases for unitary evolution over circular paths 
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Now for A£, 

tan 7g 4 + 3(l-r 2 )(A£) 2 



tan 7 g 4r 2 
For pure states (r = 1), the ratio (p4[) reduces to 1. The ratio of the visibilities, 



(37) 



v x Wcos 2 k + r 2 sin 2 k \J cos 2 K + (1 + x 2 ) esc 2 £ sin 2 K 

— — — , — — , , (38) 

v v/1 - r 2 + 2 X 2 + (r 2 -l) cosh 2tt X ^ + 2% 2 - cosh 2ir X + 2(1 + x 2 ) esc 2 £ sinh 2 ^ X 

is given in Fig. 0. Both the ratios ([34]) and (|3q) are unity for pure states (r = 1). 




FIG. 2. Ratio (g8|) of visibilities for unitary evolution over circular paths 

C. Unitary evolution of Gibbsian fc-level density matrices 

In the final part of the paper pH , Uhlmann also considered the analogous unitary circular evolution (0(3)-orbits) 
of fc-level Gibbsian states of the form 

oPrtl! 

(39) 



trace e Q n J 

where the vector J of angular momentum operators lives in fc-dimensions. (For k = 2, the model is equivalent to 
that studied in the previous section.) In Q we have conducted a detailed analysis of this model in comparison with 
the results that would be obtained implementing the methodology of Sjoqvist et al. One unanticipated outcome was 
that the Sjoqvist et al model appeared sensitive to the bosonic (odd k) or fermionic (even fc) nature of the fc-level 
systems. As an example of this phenomenon, we exhibit two additional plots (Fig. || and |]) based on analyses for 
k = 2, . . . , 11. In each of these plots, the inverse temperature parameter a is held fixed at 2. In Fig. [|, based on the 
Uhlmann methodology, below £ .65 the ten curves are all monotonically decreasing. The order of dominance among 
themselves of these curves is the natural order k — 2, . . . , 11, with the curve for fc = 2 being the dominant one. On the 
other hand, using the procedure of Sjoqvist et al, the ten corresponding curves are separated into two clusters of five 
(below £ w .44) (Fig. ||). In the first cluster (having the higher values of 7) in decreasing order of dominance the five 
curves correspond to k = 2, 4, 6, 8, 10, while in the second group (having lesser values) the order is fc = 3, 5, 7, 9, 11. 
(For other related Figures, see 0, sec. VIII].) 
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FIG. 3. Geometric phases, based on the Uhlmann methodology, for the fc-level Gibbsian density matrices, holding the inverse 
temperature parameter a fixed at 2. Below £ « .65, the curve for k = 2 dominates that for k = 3, which dominates that for 
it = 4, . . . 




-3f NN N NN\N NN 

FIG. 4. Geometric phases, based on the Sjoqvist et al methodology, for the fc-level Gibbsian density matrices, holding the 
inverse temperature parameter a fixed at 2. Below £ sa .45, the upper (even k) group of curves in decreasing order of dominance 
correspond to k = 2, 4, 6, 8, 10 and the lower (odd k) group to k — 3, 5, 7, 9, 11 



Also, we have found |13| that the (first-order) holonomy invariant in the Uhlmann approach 12, eqs. (33), (43), 
(47)], that is the trace of 



W(2ir)W(0)^ = (_l)SWp(0) 1 /a e 2iri(eoaflJ,+arinflJ w ) / ,( jl/a ) 



(40) 



reduces to (— l) fe+1 times the holonomy invariant Jl], eq. (15)] in the Sjoqvist et al analysis if we simply set the 
parameter a = sech^ in ( fl0| ) to zero, that is the influence of the angular momentum operator J y is nullified. Relatedly, 
if p(0) were the density matrix of a pure state, then the term asm9J y in the invariant ( [40| ) would drop out, with the 
trace of the resultant expression being the corresponding Berry phase |l2|, eq. (43)]. This fact was employed in |l3| ] 
to compute the Sjoqvist et al mixed state holonomy for the fc-level Gibbsian density matrices, by taking the average 
(using the eigenvalues of p(0) as weights) of these Berry phases for the pure states that diagonalize the initial density 
matrix |l], eqs. (10), (15)]. 



III. CONCLUDING REMARKS 



The analyses presented above clearly demonstrate a general inequivalence between the "experimental" approach to 
mixed state holonomy of Sjoqvist et al Q] and the "purely mathematical" scheme of Uhlmann |ll]Jl2]| . In particular, 
in the third of our analyses (sec. |ll Cf ), concerned with (9(3)-orbits of fc-level Gibbsian density matrices, it has emerged 
that the analyses of Sjoqvist et al (Fig. ||) are sensitive to the parity of k, while those of Uhlmann (Fig. ||) are not 
|p3| . (For an application of |Q to the detection of quantum entanglement, see For a further discussion of the 

results in Q, see ]l0|, sees. 2, 4].) 

A problem that remains to be fully addressed is the conceptualization of a physical apparatus that would yield 
the geometric phases and visibilities given by the analyses of Uhlmann [jll 12|. In their study jy, Sjoqvist et al did 
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describe the use of a Mach-Zehnder interferometer to "test the geometric phase for mixed states and provide the 
notion of phase and parallel transport for mixed states in a straightforward way" . However, they did not specify such 
a physical scheme for their conceptual alternative based on a purification procedure (cf. jl(| sec. 4]). 

Let us also note that the connection used for parallel transport in the Uhlmann (Bures) case is but one of an 
(uncountable) infinite set of connections ("defining reasonable parallel transports along curves of density operators"), 
which correspond in a one-to-one fashion, with the monotone metrics on the fc-level density matrices, and are all 
compatible with the purification procedure |?],[l7 18 1. In this generalized framework, the Uhlmann (canonical) connec- 
tion corresponds to the well-studied Bures or minimal monotone metric [p|,[L9|. A yet apparently formally unsettled 
question is whether or not the Sjoqvist et al connection eqs. (16)-(18)] itself corresponds to some such monotone 
metric, thus falling within this generalized scheme of Dittmann and Uhlmann of connections (and metrics) respecting 
standard purification pi. From the results presented above, it is clear that it does not correspond to the Bures metric. 
(In the Uhlmann methodology, "proper liftings" W(t) are those that minimize the length functional J Tr [dW * dW] 
of the curve of density operators undergoing unitary evolution. For a proper lifting this length becomes equal to the 
"Bures length" . One particularly interesting property of the Bures metric is that its volume element can be employed 
to determine analytically exact probabilities that a pair of quantum bits is [a priori] classically correlated |po|| .) 
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